Monomial integrals on the classical groups 



o 

(N 
o 

o 



T. Gorin and G. V. Lopez 
Departamento de Fisica, Universidad de Guadalajara 
Blvd. Marcelino Garcia Barragan y Calzada Olfmpica 
44840 Guadalajara, Jalisco, Mexico 

February 2, 2008 

, Abstract 

This paper presents a powerfuU method to integrate general monomials on the classical groups with respect 
to their invariant (Haar) measure. The method has first been applied to the orthogonal group in [J. Math. 
Phys. 43, 3342 (2002)], and is here used to obtain similar integration formulas for the unitary and the unitary 
symplectic group. The integration formulas turn out to be of similar form. They are all recursive, where 
the recursion parameter is the number of column (row) vectors from which the elements in the monomial 
^\ • are taken. This is an important difference to other integration methods. The integration formulas are easily 

implemented in a computer algebra environment, which allows to obtain analytical expressions very efficiently. 
Those expressions contain the matrix dimension as a free parameter. 

^! 1 Introduction 

I With the classical groups, we mean the orthogonal group 0{d), the unitary group U{d), and the unitary sym- 
plectic group Sp(2d) [1]. They all possess a unique invariant measure, the "Haar" measure [2], which is the 
integration measure commonly used. 

Integration formulas for the classical groups are of interest in various fields of mathematical physics. A 
number of different classes of integrals have been studied, among those: generating functions, such as the 
Harish-Chandra-Izykson-Zuber integral [3, 4], and monomial integrals, which are the concern of the present 
I work. In some cases, the corresponding integral can be solved by character expansion due to Balantekin [5] and 
■ Balantekin and Cassak [6]. Integration formulas for certain simple monomials have been developed by nuclear 
CO I physicists [7, 8]. This work was motivated by the fact that statistical methods based on the classical groups 
are very successful in describing certain aspects of nuclear reactions (see [9, 10] and references therein). Later, 
Mello and Seligman devised an algebraic method to compute low order monomial integrals on U{d) [11], and 
. Samuel solved the problem in full generality [12]. The result was an explicite formula for arbitrary monomial 
integrals over U{d). The method in [12] is based on the representation theory of U{d), and we call it the group 
theoretical method. Some 25 years ago, it has become clear that group integrals on the classical groups play an 
k> " important role in mesoscopic transport [13, 14, 15], in quantum chaos [16, 17], and more recently in some aspects 
i of quantum information and decoherence, e.g. [18, 19, 20]. 

I These applications lead to a renewed interest in efficient methods for the analytical calculation of monomial 

integrals. An interesting unconventional approach has been devised by Prosen et al. [21], and very recently, an 
improved invariant method has been developed for monomials in U{d) [22, 23], and also in 0{d) [24]. Finally, 
Collins and Sniady presented a group theoretical approach, which allows to compute monomial integrals over all 
three classical groups [25]. 

In [26], a further method has been developped, which is very different from the previous ones; the column 
vector method as it might be called. It lead to a new type of integration formula for monomial integrals on 0{d). 
This method is easily implemented in a computer algebra language. It allows to compute arbitrary monomial 
integrals very efficiently, where the result will always be a rational function in the matrix dimension of the 
group. The purpose of the present paper is to apply that method to U{d) and Sp(2d). In section 2 the matrix 
representations of the classical groups are introduced, and the notation for the different group integrals used is 
defined. The general idea of the column vector method, as well as the result for 0{d) is reviewed in section 3. 
In the sections 4 and 5 the integration formulas for the unitary and the unitary symplectic group are derived. 
A summary is given in section 6. 
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2 General considerations 



In this section, we introduce the fundamental matrix representations of the groups under consideration. Then, we 
define the integration measures and specify some notational conventions related to the integrals to be calculated. 

2.1 The orthogonal group 0{d) 

An element w £ 0{d) is a d dimensional square matrix with real entries (wij), where 1 < «,i < rf. In addition, 
w fuUfills the following orthogonality relations: 

JV 

w vF = vF w — 1 <^ Vl<i,j<(i: Wjk — 5ij . (1) 

fc=i 

The most general monomial of matrix elements of w is denoted by M.i{'w). The subscript 1 will be replaced by 
(3 = 2(4) in the unitary (unitary symplectic) case. Here, (3 is the symmetry parameter of the respective group. 

P d q 

M^{w) = \{ wi^j^ = n ' ^-^'^ = E ^^^^ ^iJ^ ' (2) 

r— 1 ijj — 1 r— 1 

where q is the order of the monomial, and / and J are multi-indeces of dimension q. The notation with the 
multi-indeccs /, J is the most common one, but in the present approach the matrix notation is more convenient. 

For the integration on 0{d), the normalized Haar measure is used. For the classical groups, in general, it will 
be denoted by ap, where /3 is the symmetry parameter introduced above. In the present case, the integral of an 
arbitrary monomial is denoted by 

(M) = / dai{w) Miiw) . (3) 



The normalization is such that (o) = 1, where o is the matrix with zero elements everywhere. 

The column vector method has been applied to the orthogonal group, first [26]. The results obtained therein 
are reviewed in section 3. It allows us to introduce the general idea of the method, as well as a number of 
notational conventions. The latter will be useful for the unitary and unitary symplectic case, also. 

2.2 The unitary group U{d) 

An element w € U{d) is a c? dimensional square matrix with complex entries (wij), where 1 < i, j < d. In 
addition, w fuUfills the following orthogonality conditions: 

d 

w = w ^ 1 Vl<i,j<(i: Wik w*i. = Sij . (4) 

k=l 

The most general monomial on the unitary group depends on matrix elements of w and . It will be denoted 
by M2{'w): 

M2iw) = n wi^j^ f[ w},j, = n 4'"' = E ^^ir Sjj^ N,, = ^ Su^ 5,j, , (5) 

r=l s— 1 ij — 1 r— 1 5=1 

where /, J and /', J' are multi-indeces of dimension p and g, respectively. Again, while the notation with the 
multi-indeces /, J and /', J' is the most common one, the matrix notation with M and N is more convenient for 
our purpose. 

For U{d), the normalized Haar measure is a2{w), while the monomial integral is denoted by 



(N\M) = J da2{w) M2{w) . (6) 
The integration formula for the unitary group is derived in section 4. 
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2.3 The unitary symplectic group Sp{2d) 

The unitary symplectic group may be defined as the subgroup of U{2d) which is invariant under the antisymmetric 
bilUnear form 

Z' = (^_°^ : Sp{2d) = {we U{2d) I w'^Z'w = Z'} . (7) 

In order to fulfill this invariance condition, the matrices w must be of the form 

' - :) . 

whith complex d-dimensional square matrices w and z. For such matrices, the unitarity conditions become: 

(w^lwi.) + (z^l^^) = 5,.,,. (i^^IzZ?^) - = , (9) 

where and i*^ denote the respective column vectors of the matrices w and z. This parametrization will be 
used to perform the integration over Sp{2d). At this point we may note that 

Vw e Sp{2d) : u)-'^ =w^ = ^^^^ t^at ^^'^ = ^"'^ = {^o 1 .) ' ^^^^ 

where w'^ is again in Sp{2d). The most general monomial on the unitary symplectic group contains matrix 
elements from four matrices: w\ z, and z^. It is denoted by: 

d 

M,{w) = M,(w, z) = n ^f' («^^)^'^ (4)"^'^ ' (11) 

where M,M' ,N,N' are d-dimensional square matrices with non- negative integer entries. If we collect these 
matrices in an appropriate way in the 2x2-matrix, we may equally well write: 

/at/ '^^ ^ 

ij — l — 1 

For Sp{2d), we denote the normalized Haar measure by (T4^{w, z). The monomial integral is denoted by 

(M) = / A<7i{w,z) Ma{w) . (13) 



3 Integration over 0{d) 

For any of the classical groups, the Haar measure is invariant under left- and right-multiplication with a fixed 
group element u, and under taking the inverse. For 0{d) these properties read: 

dai(w) /(w) = j d(Ti(w) f{uw) = j d(Ji{w) f{wu) = J dai{w) /(w"^) , (14) 

where f{w) is an arbitrary integrable function of the matrix entries of w. For a monomial integral (A/), as 
defined in equation (3), those operations translate into corresponding operations on the integer matrix M. They 
can be used to bring the monomial integral into a more convenient form. In other words, we use transposition 
and/or column permutations to collect all non-zero elements of M in the first R < d columns, minimizing R. 

The general idea is to write the integral in terms of the full d^-dimensional Euclidean space of all matrix 
elements of w, and to implement the restriction to the group manifold by appropriately chosen 5- functions. Since 
the monomial A4i{w) contains matrix elements from only the first R column vectors of ui, the integration may 
be restricted to those, as will be shown with the help of the final result, below. Meanwhile, one may as well 
assume that R ~ d. Denoting the column vectors by wi, W2, ■ ■ ■ , wr, one may write: 

(M) = /d.r(^) n-?=^|Sfv, (15) 

^1 (o) 

M''\m) = fi I / df}(«;4)n«;f/'4 n ^((^mN.)) / dn{w^) « n I / d^,,|<5(||«;dP - 1) • (16) 

f=l [■^ i=l ) u,<P •' i=l ^ 
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The matrix o in equation (15) is a d-dimensional square matrix, which contains only zeros. The subscript in the 
symbol indicates the symmetry parameter /3, which is equal to one in the orthogonal case. In equation (16) 
the integration region is the product space of R unit spheres with constant measure dr2(it;^). The orthogonality 
and the normalization is implemented with the help of appropriately chosen (5-functions. 

To convince oneself that the right hand side of equation (15) really yields the Haar measure, it is sufficient 
to realize that for R = d the measure used in equation (16) is indeed invariant under transposition as well as 
left- and right-multiplication with any other fixed group element. Similar arguments also apply in the case of 
the unitary and the unitary symplectic group. 

To evaluate the integral in equation (16), we integrate over the last column vector wn. This integration can 
be done in closed form, and the result is a linear combination of monomials in the matrix elements of the first 
R — 1 column vectors. This means, the resulting expression will be of the following form: 

AAi^«^ (Af) = ^^^^) -^i"^"'^ + ^i"^ (o) = -^i"^"'^ (o) > (17) 

K 

where K is an integer matrix with non-zero elements in the first i? — 1 columns only. The sum runs over a 
finite number of such matrices, and for [M + K) only the first R—1 columns (of M and of K) need to be 

taken into account. The second relation in equation (17) takes care of the proper normalization of the measure. 
The ratio between both equations, yields the desired recurrence relation: 

(M(^)) = ^ C{K) (M(^-i) + K) . (18) 

K 

We will eventually use the supscript of the form M^^^ to indicate that only the first R column vectors of the 
respective matrix arc taken into account. An integral, such as the one in equation (15) will be called "i?-vcctor 
integral". With only minor changes, this terminology is also applied in the unitary and unitary symplectic case. 



The one-vector integral The starting point for the recurrence relation is the one-vector integral. It has first 
been computed by UUah [8]. In our notation, it reads: 

(m) = / dai(ul) = (i),:;,ri (i),„./, - (19) 

i=l i=l 

where {z)a = r(z + a)/r{z) is the Pochhammer symbol, and rli ~ X^iLi '^i- The integral (m) is non-zero only 
if all components are even numbers. 



The J2- vector integral The recurrence relation for the i?- vector case has been obtained in [26] . It relates the 
desired _R-vector integral to a linear combination of (i?— l)-vector integrals. To obtain an explicite result, this 
recurrence relation must be continued until i? = 1 is reached, where the one- vector result can be inserted. 



{M) = Xm:e.enJ2 U ^(f 'f^f'^) E("^-«l^^) (^-) {M^"-'^ + K) , (20) 



K 



where B{a,b; z,, z,) ^ ^J^' ^^'Y' ■ (21) 

(Zi - Z2)a 

In this formula, m is the i?-th column vector of M and TTi is the sum of all its entries. The prefactor Xm:cvcn is 
one if TO is even, otherwise it is zero. The "binomial" of the integer vectors to and k is just a short hand for the 
product of binomials of their respective components: 



The symbol, (k), denotes a one-vector average, as defined in equation (19). Due to those two quantities, the 
summation is restricted to such k for which Vi : < < to^ , Ki :even. Analogous to to, k denotes the sum 
of all components of the integer vector k. The second sum in equation (20) runs over the d-dimensional integer 
matrix K with non-negative entries in the first R—1 columns, only. The following "vector-multinomial" is again 
a short hand, which reads: 

d 

{rh - k\K) = ]^(to,; - Ki\Kii, . . . ,Ki^R^i) , (23) 



=1 
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where the row sums in K must be equal to to^ — The vector fees denotes the vector of column-sums of K . It 
has at most i? — 1 non-zero components. Since our formula contains the one- vector average of fees, the summation 
over K should be restricted to those K, which have only even column sums. The last term is a (i?— l)-vector 
average, where the matrix entries of K have been added to the matrix entries in the first i? — 1 columns of M . 

As it turns out, the function B{a,b] zi, Z2) appears again in the integration formulas for U{d) and Sp(2d), 
although with slightly different arguments. 



Zero column vectors With the recurrence formula (20) at hand, it is straight forward to prove that the 
restriction of the integration in equation (16) to the first R column vectors of w is valid. This follows from the 
fact that the recurrence relation yields (A/) = (Al'---^~^^) if all components of the i?'th column vector of M are 
zero. 



Vanishing integrals It is well known {e.g. reference [10]), that the integral (Af) vanishes if any column sum 
or row sum of M gives an odd number. This can be directly seen from equation (20), which requires to to be 
even in order to yield a non-zero result. Due to the invariance properties of the Haar measure, this statement 
applies to any column- or row-vector. 



4 Integration over U{d) 

For the unitary group, we denote the normalized Haar measure by a2. It has analogous invariance properties as 
CTi of the orthogonal group, discussed in the previous section. In the present case, one finds for a fixed group 
element u E U{d): 

do'2('w) /(w, w^) = / dcr2{'w) f{uw,w^ u^) = / da2{w) f(wu,u' w^) = / d(72{'w) f{w\w) , (24) 



where /(w, z) is an analytic integrable function of the matrix entries of w and z. These operations translate into 
corresponding operations on the integer matrices M, N, which leave the monomial integral {N\M), equation (6), 
invariant. In particular: (i) simultanous column permutations: {N\A4) = (7V7r|A'/7r), (ii) simultaneous row 
permutations: {N\M) = (ttNIttM), where tt is an arbitrary permutation matrix, and (iii) conjugate transposition: 
{N\M) = {M'^\N'^). Due to the invariance under column permutations, we may assume without loss of generality 
that the non-zero elements of M and A'' are all restricted to the first R < d columns. In other words, this means 
that the monomial A^2(w) contains no matrix elements from column vectors Wp, with /i > i?. 

Analogous to the case of the orthogonal group, we again write the integral in terms of the flat Euclidean 
space of all complex matrix elements oi w £ U{d). Then, we implement the restriction to the group manifold by 
appropriately chosen ^-functions: 

(AIM) = / da2H n ««)"'^ - ^^^W^ ' (25) 

J ^^tl Mr {0,0) 

Af^''\M,N) = n ( / df72(«j^)n<-^ i^^rA n ■ (26) 

The subscript 2 in the symbol J^2^'' stands for the unitary case (3 — 2. The first R column vectors of the unitary 
matrix w are denoted by wi, . . .wr. The integration region in equation (26) is the product space of R unit 
spheres with constant measure di}2iwi^): 

J dn2(.w0'^jiy d'iw,O^Si\\w^f-l), (27) 

where we define the flat measure on the complex plane via z g C : d'^{z) = d(Re z) d(Im z). The (5-functions 
in the equations (26) and (27) implement the orthogonality conditions and the normalization. Note that the 
(5-function in equation (26) really is the product of two (5-functions: S^'^\z) = (5(Rez) 5(Imz). 

Since we assumed that the monomial M2{w) contains matrix elements from only the first R column vectors, 
we may restrict the integration in (26) to those vectors. As in the orthogonal case, this does not affect the result 
of the monomial integral, as will be shown below with the help of the final result. Meanwhile, one may as well 
assume that R = d. Note that for R — d the integration measure in equation (26) is indeed invariant under the 
transformations performed in (24). This guarantees that equation (52) really yields the Haar measure. 
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4.1 The one- vector formula 

In the one-vector case, the matrices M and N can be replaced by the d-dimensional vectors m and n. In that 
case, there are no orthogonahty conditions to obey. We may write: 

(n|m) = / da,{w) {{w^^' ^^^J,^ Ni'\rn,n) = / d.n,{w) {{w^ . (28) 

i=l ^^2 [0,0) J 

As suggested in equation (27), we integrate over the full space R^'', while the normalization is implemented with 
the help of a 5-function. This introduces an integration constant, denoted by C2{d, 1): 



M^'\rn,n)^C2{d,l)f[y d\w.) «)"' | <5 (HuJf - l) 



(29) 



The J- function is removed as follows: Setting Wi — Ui/^ we get: 



A/-W(m,n)r'^+(*+«)/2-i = C2(d,l)n| I d'^uT' «)"' | <5 (E. " ' (30) 

where fh ~ ^« ^^^^ ^ ^ ^i^i Multiplying both sides of the equation with e~^ and integrating on r 

from to 00, the (5- function disappears: 

d ^ d 

U^'\rh,n)r{d+^^) =C2(d,l)n / ^^M^r «)"■ c-l"-'' =[]/(m„n„0) . (31) 

The Gaussian integral /(m, n, a) is defined in equation (92) in the appendix. The general result for that integral 
is given in equation (100). For a = 0, it evaluates to 

/(m, n,0) = TT (5„,„ n! . (32) 

From this, it follows: 



M^^\rn,n)T{d+^^) = C2(d, 1) H A^f ^ (o, 5) r(d) = C2(d, 1) tt^ (33) 



i=l 

where (5m. n denotes the product of Kronecker deltas between the components of ffi and fi. This finally leads to 

d 



{d)-^ X{m,\ ■ (34) 



4.2 The i?- vector formula 

To solve the i?-vcctor integral, we start from equation (26) and separate the integration over the first i? — 1 
column vectors from the integration over the last clumn vector wji. Let us denote the i?-th column vectors of 
M, N and w with m, 72, and respectively. 

H-i r » d ^ R-l 



Ni''\M,N) = n / df^2(^e) t{wf' [wl^f^^ n {-^^'^((^mN.))} 4''\m.n) (35) 
5=1 {.•' i=i J ti<v 

J-f )(m,n) = /■ df72(^)n{<- «)"■} n ■ (36) 

1=1 /i=i 

We start by flattening the integration measure f22, using the technique from the one- vector integral. With 
d^(u)i) ~ d{Kcwi) d(Imz«i) we may write: 

J^f )(m,n) = C2{d,R) n I I d\w,) «)"' | ^ (E. ^ 1) H -^''^ (Hl^)) " (37) 

1=1 '^•^ ^ u=l 



6 



The transformation Wi = Uij^/r leads to: 



11 ^ ^ i?,— 1 

1=1 ^ ^l=l 

Mulitplying both sides with e^^ and integrating on r from to cx3 gives: 

J^''\m,n)r{d-R+1 + Si±^) ^C2id,R) J] / d^^) K)"' e-<"4 J^^^ ((idjw)) . (39) 

Using the Fourier representation of the one-dimensional delta function, we may write: 

(40) 



With the help of this representation, we write: 

R-i fl-i . ,2 



n ^^'^ (s. "0 = n / ^ exp {r, - 

Therefore 

R-l r ^ ,2/ N ^ d -R-i 



(41) 



V 2 / /i^l ^ i—l /J,— 1 

/(m,7i,a)= /"d2(u) u™ (u*)" e- 

^.-rt,-.. (::„:,. ..... 



v=0 



For details about the computation of that integral, see appendix. In order to put the product of the functions 
f(mi, rii, ai) into a suitable form, we define additional integer vectors k and p, where pi = min(TOi, n^): 



d P d f 

n/(™.,n.,«o=^' {-r E(-)'n{ 

i—l K.—0 i—l 



K,\ (a*)"-""- a^""' \ e-2:j".|- . (43) 



With m and k, we denote the sum of vector components in rh and k, respectively. For the expansion of the 
powers of the coefficients ai and a*, we need the integer matrices K and L. Both are c?-dimensional matrices 
with non-zero elements in the first R — \ columns, only. 



d 

rrii—Ki . ^ _ 



i=l i—l /C /J.— 1 2—1 

_R-1 d 

= ^(m - ^|X) n (^;)'" n ' (44) 

where fci, . . . fcij-i denote the column vectors of K , while the overbars denote, as usual, the summation over all 
vector components. In a completely analogous way, we write also: 

d R-l d 

n^"""* = - ^1^) n n^)"^"" - (45) 

z— 1 L fi—1 i—l 

where /i, . . . Ir-i denote the cohimn vectors of L. Thus, we obtain from equation (42): 
J'f)(77i,n) 



r(fl)/-=. C2(d,i?)^'^(-)'^ ^fiA^^^A^" 



r(d-R+i+si±^) T^y^^ n I 

V Z / K, — l—l 

J2im~K\K) (n-«|L) / i^'P"" c'^.l"'!^ [J^;,^- «)^- . (46) 

rv' r , 1 I' 1 



K,L ti = l • 
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Due to the fact that the integration in equation (26) is restricted to column vectors zZJ^ which are pairwise 
orthogonal, we obtain: 

51 = 5Z 5Z '^^ XI ^^'^ = H I^A"!^ ■ (47) 

This is a crucial point in our calculation. It eliminates the remaining matrix elements of w from the exponential, 
and J'2^\rn,n) turns into a linear combination of simple monomials. The integrals over are now easily 
evaluated, using equation (31) from the calculation of the one- vector average: 



?{r^::^; B-''n{(:)(::)^ 



R-l 



xJ2{m- k\K) {n k\L) [] 6-,^j^ k,\ J] W^)'"'' ■ (48) 

K,L ju=l i=l 



Now, we insert this result into equation (35) and write for {M, N) and M!^^'^ (o, o) : 

M=l £ = 1 I 2=1 J Im<'^ 



mi\ in. 



^(m-^|/0 (n-^|L) n {^fe„r, fc^i} AAf"i)(Af(«-i)+X,iV(«-i)+L) (49) 



R-l 

X ' - ^' -'^ 

K.L /i=l 



(fl) C2(rf,i?)^'^-^+l ..iR-l). , 



Thus we obtain: 



{N\M) = (-)- id-R+ 1)^^^,2 ±i-r n { (I:) 

x^(m-K|A-)(n-.J|L) J] {^^.^^a^'} (^'''-'^ + + A') • (51) 



Note that due to the multinomial coefRcients (m— k\K) and (n— k|L), for the column sums of K and L, it holds 
that: yZ^Zl fc^ = m — K, and also X]^=i l^j, = n ^ k. In addition, to obtain a non-zero contribution, the column 
sums of K and L must be the same also: Vl < fi < R — 1 : ~ l^. Together, both conditions imply that for 
fh ^ n, the whole monomial integral gives zero. Therefore: 

' m\ / n 



K — O 

X 

K,L 



{N\M) = Srn.n E (^j y (^l'^) i?(m,«;d,i?- 1) 

^(?7i - k\K) {n - n\L) {Qhs) (iV^^^i) + L\M'^^-^^ + K) . (52) 



In this formula, we use the vector notation for the product of binomials and multinomials, as introduced in 
section 3. We also express the products Y\^^i Kil and n^=i i^, ^tJ^^- terms of the one-vector averages, given 
in equation (34). In this way, the coefhcient B{fh,R;d, R — 1), defined in equation (21), comes into play. The 
summation over k is restricted to those vectors, for which Vi : < < min(m.i, rii). K and L are c?-dimensional 
matrices with non-negative integer entries, which have non-zero elements in their first R — 1 columns, only. The 
summation over K is restricted to those K, for which Vi : Ki^, = mi — Ki. Similarly, the summation over L is 

restricted to those L, for which Vi : L^^ = rii — Ki. In addition, the one- vector average (^csl^cs) implies that 
the column-sums of K must agree with the corresponding column sums of L. These column sums arc arranged 
in the first R — 1 components of the d-dimensional vectors kcs and ^cs, respectively. 



8 



Zero column vectors Assume that in the monomial integral {N\M), the i?'th column of M and N contain 
both only zeros. Equation (52) then yields: {N\M) = (iV(^-i)|M(^-i)). This proves by a simple induction 
argument, that the integration space in (26) may indeed be restricted to those column vectors, which contain 
matrix elements from the monomial M2{w). 

Vanishing integrals Denoting monomials in the form: wJl^j,^ • ■ • '"'I b ''^aipi ■ ■ ■ Wa^p^ (see for instance, refer- 
ence [14]), it has been shown that the monomial integral (it;* . . . w* ^ Waifii ■ ■ ■ Wa^fj^) vanishes, unless p = q 
and ai, . . . , ttg is a permutation of ai, . . . , Op and /3i, . . . , /3g is a permutation of &i, . . . , 6p. This is equivalent to 
the requirement that all columns sums and row sums of M agree with the corresponding column sums and row 
sums of N. In view of the invariance of the Haar measure under row- or column permutations, this is a trivial 
consequence of the leading factor (5m, n in equation (52). 



5 Integration over Sp{2d) 

For the unitary symplectic group, we denote the normalized Haar measure by a^. For 



~w z 



Sp{2d)3w^[ \ , let .f{w)= f{w,z,w*,z*) , (53) 



where /(wi, zi, u'2, Z2) is an arbitrary analytic(?) and integrable function of the matrix elements of wi, zi, u'2, 
and Z2. The invariance properties of the Haar measure then imply: 

da^^w, z) f{w) ^ / dai{w, z) f{wu) = / dai{w, z) f{uw) = / dcri{w, z) f{w^) , (54) 



where it is a arbitary but fixed element of Sp(2d). 

Since we are interested in the case, where the function / is a monomial, we use some particular permutation 
matrices from Sp(2d) to obtain invariance relations for the monomial integral (M), under certain transformations 
of M. For the following considerations, let 

w = ( ^ ^ ] , u = ( ^ ^ ] . (55) 

(1) For It = TT, a (f-dimensional permutation matrix, and v = o, multiplication from the left leads to 

/ — IV 7T 7r \ 

w ^ wu ^ [ * , (56) 



-ZTT —W IT 



such that 

d 



where M' and N' are the sum of all matrix entries in AI and iV, respectively. Similarly, multiplication 
with u from the right yields: 

(M) ^ / d.. n/«)«"(-"")f'(«*)S"(-''«'*)?' - ((:-".«. ;-"m)) ■ (^«) 

(2) For u ~ o and v ~ tt, multiplication from the left yields: 

m ^ / d.. n/-)"^^ (-)f = ( "X-^^ m 

whereas multiplication from the right yields: 
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(3) Finally, replacing w with its inverse, w\ yields: 

(M) - Jda,{w) l[{-w,.rH^^hrH-w^^^^^^^ (61) 

where M and TV are the sums of all matrix entries of M and N, respectively. 

As we will show below [sec equation (88)], the monomial integral (M) is zero, unless M + M' = N + TV'. 
Equation (61) also implies that (M) is zero, unless M + N' = N + M' . Both conditions together imply that 

(M) 7^ ^ M = M' = N ^ N' . (62) 

We can therefore savely ignore the sign-prefactors in (57), (58), and (61). The results in item (2) prove the 
invariance of {M) under synchronous column permutations [equation (59)] and synchronous row permutations 
[equation (60)], applied to all submatriccs of Af. Therefore, the equations (57) and (58) add only two additional 
operations: 

'N' M\ (M N'\ fN' M\ fN A/'^ 

M'l ^ [m' N j ' [ N M'j ~^ [n' M 

The integration 



(63) 



It is again possible to write the integral over Sp{2d) in terms of an integral over the flat Euclidean space of all 
complex matrix elements of w and z, implementing the restriction to the group manifold by appropriately chosen 
(5-functions: 

(M) = I da,{w, z) M,M = ^^^5^ , M,{w) = ^ 4'"' ««)'^*' ■ (^4) 

(o) 1,4 = 1 

Here, we assume that the monomial M.ii{w) contains no matrix elements from the column vectors w^,z^ with 
fi > R. One may have used the invariance of (M) under the above discussed transformations in order to minimize 
R. The matrix o is of the same shape as M, but all its entries are equal to zero. 

AAf )(M) = I n / d"4(^/??, %) i M,{w) n S'^'Hi^.l^u) + {z,\z.)) S'^'Hi^,]^.) - HWr^)) • (65) 

The subscript in the symbol is the symmetry parameter /? = 4, since we consider here the unitary symplectic 
case. The first R column vectors of the matrices w, z are denoted by zwi, . . . w^, and zi, . . . zr, respectively. The 
integration region in equation (65) is the product space of R unit hypersphcrcs with constant measure drJ4(i(;^, z^): 



y" df}4('(«5,z^) oc d'^(wi5,Zj^)| (JdltZljl 



|r + ||z^||'-l) w,zeC : d''(«;,z)=d^(«;)d^(z), (66) 

where d^(z) denotes the flat measure on the complex plane, as already used in equation (27). The (5-functions 
in (65) and (66) implement the orthogonality conditions and the normalization, as defined in equation (9). 

According to the assumption, the monomial M,^(w) contains matrix elements from the first R columns of 
vj and z, only. Therefore, we restrict the integration in equation (65) to those column vectors. The fact that 
this does not affect the result of the integral, will be shown only at the end of this section, where we have the 
final result, equation (89) at our disposal. Note that the integration measure in (65) is indeed invariant under 
the transformations performed in (54). This guarantees that the above construction really implements the Haar 
measure. 



5.1 The one- vector formula 

In the one- vector case, the monomial M.t^(w) contains matrix elements from one column of w only. This means 
that the matrices M, M',7V',A^ have only zero entries in all columns, except for one which must be the same 
for all matrices. We denote the column vectors with non-zero entries by rfi^rfi' and n' and use the following 
notation for the monomial integral: 



(M) 



da,{w,z) \{wrzT'^{w:r{z:r 



(67) 
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In order to evaluate that integral, we observe that there are no orthogonality relations to obey. The integral 
must therefore be equal to the one- vector integral {q\'P) over the unitary group U{2d), where p and q are the 
concatenations of the vectors m, m' and ft, n' , respectively. Therefore: 

d 

= SmM Sm',n' {'^d,)Z^\rh' H ' ^^^^ 

i=l 

where ifi and ifi' denote the sum of all components of rh and m', respectively. 




5.2 The vector formula 

We aim at a recursion formula, which relates an _R- vector integral to a linear combination of simpler (_R— l)-vcctor 
integrals. To derive such a relation, we separate the integration over the i?-th column vectors of w and z from 
the rest. Let us denote the i?-th column vectors of M, Af, A^, N' by to, to', n, n' . Then we may write: 



R-l 

J] [5^^\{w,,\w,) + {z^\z,)) S^'H{z;\w,) - {w;\z,))} J^''\m,m',n,n') (69) 



R-l 

X 



:7f)(TO,TO',n,n')= / dn,{w,z) l[{wrzf'{w*r{z:)<} l[S^^\{w,\w) + {z,\z))S^^\{z;\w)~{w;\z)). 

(70) 

We start by flattening the integration measure $^4. With equation (66) and d^(z) = d(Rc z) d(Imz), we may 
write: 



jf)(...) = c4(d,i?) njl d2(^.) 1 d'iz,)wrzfiz:)<iw*r^ 



n -Jf^HKN) + i^M) s^'H{^,\w) - ■ (71) 

The delta function for the normalisation can be removed, using similar steps as in the case of the unitary group. 
This yields: 

■ •) - r(2(.-i..l)?!^?^'..' + .)/2) n { / d^K) / 
X n + {zM) S^'H{^,\u) - {w;\v)) . (72) 

According to equation (40), the Fourier representations for the remaining delta functions may be written as 

n^^^^ + z:^v,)] = n / ^ ne^'i-[-^(-.>'+<.-)i ^ n | / ne^'i™(--+^-) 

fj,= l /j=l i=l /j=l '-•^ ^ i=l 



fl-l if-l 

1 /I— 1 //—I /J,— 1 
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Insertion into equation (72), and the exchange of the order of integration gives: 
T{2{d- R+l) + {Th + fh' + n' + fi)/2) i-^ \J 7r2 J tt^ 



X 

1=1 
Ci{d, R) 



u=l '-•^ "^-^1=1 



(74) 



The function / is the same which appeared in connection with the i?-vector integral for the unitary group, 
equation (42), and which has been computed in the appendix. Inserting the resuh, equation (100) and expanding 
the powers of [at + 7i) and {[3i — £,), we obtain: 



p / \ / \ m — K n~ K 



K,—0 /c— Z— 

(75) 

/(^^n^/5-s) = .cH^-^|^E(-)™'-«Y3)n.'! E 5: 



V 

K'=u r=u 



where we have for the moment suppressed the index i. In order to completely expand the product of the functions 
f{mi,ni, Ui + 7i) and f(m[,n[, Pi — e^), we define additional vector indeces: pi = mm{mi,ni),p^ = min(m^, n^), 
Ki, and k'^. The complete expansion then reads: 

d 

n f{m^, n„ a, + 7«) /(m^, n^, /3, - e,) = tt^'^ e" 2:.(l".+7.P+lft-e.P) 

K — O \ / \ / f^'—g 

rh—Kfi—K. / _, _,\ /_. _.\ rn' —k' n —K 



k' J\ I 



X iiyoii) Kli) (^l li (Pj'(-ei) • ' • A (-eO * * (77) 

The argument of the exponential function can be simplified, as follows: 

X] + 7*1^ + E 1-^' ^ ^ E E [(^m"^*/. + Tf,z,^){a*w,^ + T*zl) + {af,z*^ - T^w,^,){alzi^ - t*w*^)\ 

i i f-LV i 

= E + {Zf,\z^)) + Tf,al{{z*^\w^) - {w*^\z^)) + af,T*{{w^\z;^) - {z^\wl)) + ^i^Kiiz'^K) + i^J^lwl))] 

= Y.^\o^? + \r,?). (78) 

where we have used the orthogonality relations in (9). For the expansion of the powers of coefficients a*, 7*, 
etc., we will need eight matrix indices, each of dimension dxd. It will be convenient to combine them into two 
2x2 block-matrices, as follows: 
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Then we combine pairs of coefficients, which lead to monomials of entries from the same submatrix of w: 

{-e^^-«^{r-f-^-'' m,){n'-K'-k'\K2) f[Kf H^""^'''' (80) 

i=l Ki,K2 ^=1 s=l 

t[iP:f^ ir"'-'- = E (^'1^0 (" - - UKf^' ^ n 4"^^''"'^ (81) 

d _R-1 ^ d 

i— 1 L1.L2 i—1 

n(7*)'— /3f' ^ {l'\L[)irn-^-k\L',) H aj' (r;)^.' n(z,*,)^->+^- , (83) 
t=i ^'1,^2 '^"-'^ 

where we eventually convert the subscripts 1, 2 of the eight matrix indices into superscripts, especially when we 
refer to particular matrix elements. The parameters fc^, fc^', etc., denote the column sums of the matrix indices 
Ki,K[, etc. . Plugging the complete expansion into (74), we obtain: 



(^. ^^ n, ^ ^ ^^^^ ^ I 1 1^^^ sgj'^ (J y u) u^^'-"^-^ 

rfi — K rh' — k' n — K n' — k' 

E E EE( 



X 

k—d k'—o l—o I'—o 



K—o k'—o i—1 

k'+T' I ^ \ I fTi' ^ / n — k\ / n' — k' 



V k' J \ I / V I 



X Y ik\K^)ik'\K[)if\Li){l'\L[) {rt -k' -f'\K2){n-K-r\K!2){m' ^k' -k'\L2){m-K-k\L'2) 

KuK2 



X 

^ttA Kl,+Kl K}:,+Kf„ , , 



n n v^ *" «)"-+"- (^^)"-+"'^ . (84) 

z— 1 

We can now merge each sum and its corresponding binomial from the second hne with two sums and their 
corresponding multinomials from the third line: 

Y.\%) E (^i^')("^-^-^i^2)- E ("^-^1^2,^1) 



E E (^'i^i) ("^' - - ^'1^2) = E (^^' - ^^'1^2, i^o 



n — K 

' n — K 



l=. 



E( I ) E (^i^i)(^-'^-te= E ("-'^i^i'^^) 



E C"';/') E (^^l^i)("'-^'-^^l^2)= E (^'-^'IM'^2), (85) 

where an extended vector- multinomial such as (mjifi, L'^ means that the rows of K\ and L'^ should be concate- 
nated, to yield a product of multinomials, each with 2d elements. Evaluating the final integrals on cr^ and t^, 
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we obtain: 



jl^\rn, rn' , ft, n') 



~i-\~m' -\'n-\~n' \ / j / j ^ ' 



r(2(Ar-i?+ 1) 



K^J< 
d 

X n w. 



/ 771 — 


K 


L'2 




f ft — 


K 






I m' — 


k' 


L2 






k' 


L[ 


K2J 



K. — K—O 



W,^' 2,;^^ K^) '^^ (V) 



(86) 



1=1 



where, we have replaced fc' and /' with A'( and Z'j^, as a more appropriate notation for the sum over all matrix 
entries in K'^ and i'^, respectively. We have also further expanded the vector multinomials, concatenating the 
matrices and vectors in the vertical direction. 

Now, we may insert this result into equation (69) in order to obtain A/'J^'' ( Ai") as a linear combination of 
terms which contain the (i?-l)-vector integrals A/'j^"^'' (Af'^''^"^) + 1^1 + 7^2) ■ Note that the matrix indices in (86) 
are chosen in order to fit into that scheme. Dividing through the corresponding expressions for M\ (o) and 
we obtain: 



(M) = 



r (2(d - i? + 1)) (-)™+™' 

r (2((i - i? + 1) + (to + m' + n + n') /2) 

E 



E(- 



, (2d)s+s' 



X 



The vectors fc^, fc^„ 



L'2 i^i 
L2 



n — K 
n' -k' 



L[ K2 



n 



^cs 



m — K 

X {M^^-^^ + Ki + K2) . (87) 

1 1.2 n. g^jjj^ d-dimensional vectors which contain in their first R — 1 components the 

column sums of Ki and K2 in the following order: k^^ contains the column sums of Ki and K'^] fc^g contains 
those of K2 and K2, 1^^ those of ii and L[; and 1^^ those of L2 and i2- Note that due to the one vector integral 
in these column sums, one may conclude that 



R-l 



E + = '1 + '1' + '"^ + 'I 

This allows us to use again the function B{a, b; zi, Z2) defined in equation (21): 

^ B (fh + fh' , R, + R'; 2d, 2(i? - 1)) 



X E 



/ TO 


— K 


L'2 






-k' 


L2 





n — K 
n' ^k' 



Li 

L[ K2 



{-) 



7i 



n 



{M'-"-'^ +K1+K2) . (89) 



The first line of this expression agrees precisely with the corresponding part of the recurrence relation for the 
unitary group U{2d). The differences appear in the second line, where we have to sum over 2x2-block matrices 
Ki and K2, with the restriction that the row-sums of the block-matrices in the following vector- multinomials 
must agree with the corresponding component of the given 2(i-dimensional vector. 

Zero column vectors In the case that the _R'th column vector of all four matrices M, M' , N, and N' contain 
only zeros, the summation over k and k' contains only one term, namely k ~ k' = o. Similarly, the summation 
over Ki, K2 also contains only one term, where all matrix entries of Ki and K2 are zero. Therefore, wc find 



(M) 



N' 


M\ 




\N 


M' 


}-{ 



= (M 



(B-l) 



(90) 



As in the case of 0{d) and U{d) in the previous sections, this justifies equation (65) where we have ignored 
column vectors of w and z whose entries do not occur in ^^{w). 
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Vanishing integrals The monomial integral (TVf) obviously vanishes if m + m' is different from n + n'. Taking 
into account the transformations from (54), which do not change the integral, this means that the integral 
vanishes unless the column sums of M fulfill the condition: 

V 1 < j < (i : rhj = rhd+j , (91) 

where ffij is the column sum of column j of M . An exactly analogous condition holds for the row sums. 



6 Conclusions 

We applied the column vector method^ originally developed for the orthogonal group in [26] , to monomial integrals 
over the unitary and the unitary symplectic group. As in the orthogonal case, we obtained recursive integration 
formulas, where the recursion parameter is the number of non-empty column vectors (these are those matrix 
columns of the group element, which contain at least one matrix element from the monomial). In all three 
cases, the recursion formulas are of a very similar form and relatively easy to implement in a computer algebra 
environment. This provides an efficient way to compute such integrals analytically as a function of the matrix 
dimension. 

The approach presented here is very different from the group theoretical approach developped in [25] , which 
also provides explicite formulas for the analytical calculation of arbitrary monomial integrals over the classical 
groups. We did not yet compare the efficiency of both methods, but we would expect that the result of such a 
comparison will depend on the problem at hand. 



A Gaussian integral on the complex plane 

Here, we compute the Gaussian integral from equation (42). Writing for the complex integration variable 
u — X ~\-iy we get: 

f{m,n,a) = // dxdy li™ (u*)" e""*" e" . (92) 



First note that /(m,n, a)* — f{n,m, —a), which implies that it is sufficient to consider the case m > n. Then, 
we use the complex differential operators dz and dz to write: 



/(to, n, «) = (-)" a™ as yy dxdy e-" "e""-" " 2dz^d,-idy 2dz=d^+idy. (93) 
With a — ai -\- ia2, it follows 

/(to, n, a) = (-)" a™ I dx e-(^'-2ia2 x) f ^-{y'-2ia^ y) 



= (-)" a™ a;j TT e"" " = tt (-)"+" dl (a*)™ e-" " = TT (_)'"+" P(m, n) e"" " , (94) 

where Pim^ n) is a polynomial of total order m + n, with the property that the difference between the power of 
the conjugated variable and that of the non-conjugated one is m — n, fixed. Hence, for P(m, n) we consider the 
following Ansatz: 

n n 

P{m, n)^J2 Cfc"'"' («*)™-"+'^ a'' = (a*)™-" ^ C^™'") (a* a)'' . (95) 

k=0 k=Q 

For n = 0; the polynomials P(to,0) are known, while for m > n > 0, they can be computed by recursion: 

n = : P(to, 0) = 1 => C^"'"' = 1 
TO > : P(to, n) = e"' " B'^ (a*)" e""* " 



P(to, n + 1) e"' " 82 P(to, n) e""' " = -a P(m, n) + J2 ^1™'"' {m~n + k) (a*)™-("+i)+'= a'' 

(96) 



Equating the coefficients for corresponding powers, we obtain for m < n < k < 

'{m-n) C^™'") : fc = 



l<k<n = 1 . (97) 

fc = 71 + 1 
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Ordered according to n and k, these coefficients can be arranged in a two-dimensional scheme, forming a Pascal 
triangle. From this scheme it is easily seen that 



C^"'"^ =: (m - n + fc + 1) (m - n + fc + 2) . . . m ( ) . (98) 



Therefore, wo finally obtain 



k „ n — k 



^ ' ' {m - n + k)l kl {n - ky. ' ' (m - k)] k\ {n - ky. ' 

(99) 

As mentioned before, here we have assumed that m > n. For m < n hy contrast, we use: f{m, n, a) = 
f*{n^m, —a) so that in summary, we may write: 

/(m, n,a)=T: (-)" X!^")'' (T) (fc) e""*" P = mm(m, n) . (100) 
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